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Abstract 

We present a new Ginzburg-Landau theory for superconductivity in UPts, 
based upon a multicomponent order parameter transforming under an irre- 



t> ■ ducible space group representation; the phase is staggered in real space. Our 
(N 

■ model can explain the H — T — P phase diagram including the tetracritical 

00 

■ point for all field directions. We motivate this unconventional superconduct- 

^ : 

0^ . ing state in terms of odd-in-time-reversal pairing that may arise in one- or 



two-channel Kondo models, and suggest experimental tests. 



C ■ PACS numbers: 74.20.De, 74.70. Tx, 74.25.Dw, 64.60.Kw 

o ■ 

I Heavy fermion superconductors have continued to attract attention since their discovery 

I in 1979 |T|]. In these materials, fermionic excitations with effective masses hundreds of times 
the free electron mass undergo a pairing transition. Anomalous power laws in low temper- 
ature thermodynamic and transport properties together with complex phase diagrams 
suggest interpretation in terms of Cooper pairs with exotic symmetries. For example in 
the magnetic field-temperature H ~ T plane for ambient pressure P \ bar, the hexag- 
onal material UPta apparently possesses three distinct superconducting phases, shown as 
A,B,C in Fig. 10. In the P — T plane for zero H there are similarly three observed 
phases . Intriguingly, the B phase in UPts appears to break time reversal symmetry T in 
that muon spin rotation spectroscopy reveals the development of tiny additional magnetic 
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moments (order 0.001//^) @]- So far, no comprehensive microscopic theory exists for these 
materials and so much theoretical effort has been devoted to phenomenological approaches. 
No such approach to date has successfully explained UPts without substantial fine tuning 
of parameters. 

In this Letter, we propose a new model for the phase diagram of UPta in which finite 
center of mass momentum (FCM) pairs give rise to an order parameter which transforms 
as a triplet under the operations of the full hexagonal space group. Our model produces a 
plausible description of the H — P — T superconducting phase diagram. In particular, we can 
explain the tetracritical point found for all field orientations. While the Ginzburg-Landau 
free energy contains numerous parameters, experiment constrains the zero field values to a 
narrow region which allows for a T-breaking B phase. The crystalline structure of UPta 
plays a crucial role in our considerations. The proposed superconducting order parameters 
are microscopically motivated in terms of odd-in-T pairing such as may arise in the 
two-channel Kondo model |^,|^ suggested for heavy fermion materials |p. Such pairs may 
generically give rise to negative pair hopping energy favoring staggered superconducting 
states 0. We shall discuss data supporting this odd-in-T pairing and point out crucial tests 
of our model as well as some of its limitations. We note that our theory is applicable only 
near Tc and ifc2, and is macroscopic in character. Explanation of the low temperature power 
laws awaits the development of a suitable microscopic approach. 



To place our work in context, we first discuss the E doublet model [T^. In this picture, 
the Cooper pairs have zero center of mass momentum (ZCM), and are described by basis 
functions (f)a{k), <f)b{k) which form a closed vector space under the symmetry operations of the 
hexagonal space group. With each basis function is associated an order parameter amplitude 
rja^b, defined through the gap function A{k) ~ ria4>a{k) + //^^^(fc). The phenomenological 
approach consists of writing down and minimizing a Ginzburg-Landau free energy which is 
an expansion in the f] amplitudes. A zero field splitting is present which most likely originates 
in the orthorhombic strain induced by magnetic order. Then the three phases have rja 7^ 0, 
rjb = O(A-phase), rja 7^ 0, 77?, 7^ O(-B-phase), rja = 0, rjb O(C-phase). The apparent 



time reversal breaking in the S-phase is explained by suitable a choice of Ginzburg-Landau 
parameters which energetically favor a relative phase factor of i between rja, rji,. The E model 
has a flaw: while the data show a tetracritical point for all field orientations, a tetracritical 



point generically arises only for in-plane magnetic fields [|ri|]. This problem is remedied by 
fine tuning to zero the gradient term in the free energy which mixes the a, b components 
|T0| . As a possible alternative, the authors of Ref. propose a model of accidentally 



nearly degenerate ZCM singlet states (DS model) which have no such gradient mixing term. 
While this can explain the tetracritical point for all field orientations, there is no compelling 
explanation for the accidental near degeneracy. 

Our work proceeds via a different phenomenological route: by considering Cooper pairs 
with FCM, we are free to examine space group representations which merge the desirable 
properties of the E doublet and DS models by having order parameter degeneracy and the 
absence of gradient mixing terms enforced by symmetry. 

Before turning to the core of the paper, we shall briefly review the motivation for this 
work in terms of odd-in-T pairing theory. For an odd in T pair wave function it may be 
shown within a quasiparticle framework |^,^] that (i) for ZCM pairs, the "Meissner stiffness" 
which goes as 1/Xl(T)^, \l the London penetration depth, is negative] (ii) an immediate 
corollary is that the pair transfer energy between two odd-in-T slabs, or Josephson coupling, 
is negative-alternatively, the coefficients of the gradient terms in the free energy are negative. 
To the extent that these properties are generically true (it is unknown whether a quasiparticle 
picture is applicable in all cases), we then expect a non-uniform superconducting state in 
which the pairs have FCM and the phase is modulated in real space. Indeed, the Majorana 
Fermion treatment of the ordinary Kondo lattice model favors just such a state, with the 
superconducting order parameter fully staggered, i.e., the pair hopping from site to site is 
negative 0. 

There is so far no numerical evidence supporting this odd-in-T picture for the single 
channel Kondo lattice model [|T^]. There are strong reasons to consider this state in the two- 
channel Kondo lattice model. In this model, two identical "channels" of spin 1/2 conduction 



electrons interact antiferromagnetically with a lattice of spin 1/2 local moments. For this 
model it has been shown that (i) in the impurity limit, the pair field susceptibility for 
the odd-in-T order parameter which is a singlet in spin and channel indices diverges as 
— ln(T), T 0, and (ii) all heavy fermion superconductors have appropriate symmetry and 
dynamical conditions to be describable by two-channel Kondo models [P|,P^ . For the uranium 



based materials a two-channel quadrupolar Kondo model can apply, with spin corresponding 
to local quadrupolar indices and the channel labels corresponding to conduction magnetic 
indices. The spin and channel singlet states transform as the A2 representation for each 
point group relevant to the heavy fermion materials [|14| {A2 ~ (x^ — 3xy'^){y'^ — Syx"^) for 
the hexagonal point group). 

Assuming that such local odd-in-T pairs of A2 point symmetry form at each site with 
negative pair hopping in the heavy fermion superconductors, we can determine the allowed 
multidimensional space group representations which describe the ordering in each compound. 
Note that such a description is reasonable only to the extent that the Kondo screening 
length is strongly reduced in the lattice relative to the impurity, which is in fact expected 
in the single channel Kondo lattice model |TH]. In a hexagonal crystal with one 4f/5f atom 



per unit cell the pair energy minima will be at the two K points of Fig. 2. For UPta, the 
crystal structure with two atoms per unit cell and interlayer frustration for antiferromagnetic 
(phase) couplings pushes the pair hopping minima away from the K points, and for a wide 



range of parameter values they move to the M points which have a three point star |T6|. In 
addition, the presence of two atoms per unit cell gives rise to non-degenerate bonding and 
non-bonding combinations of pair orbitals within each cell. We assume that the non-bonding 
combination is sufficiently high in energy to be negligible in the current analysis. 

We now turn to the discussion of a specific Ginzburg-Landau model for UPts based 
upon the 3 dimensional M-point representation. While motivated by the odd-in-T order 
parameter discussion above, we are not limited to this as a potential microscopic source for 
our phenomenology. Free energy invariants are constructed by standard group theoretical 
methods, and we have checked the completeness of our set using image group techniques [jl7 



We include a linear coupling to an orthorhombic strain field e, which breaks the hexagonal 
symmetry. The free energy density is given by 



F = Fo + Ff + (la) 
Fo = a+|r/3r + "-(l^iP + \V2?) + A E 

n 

Pf = Y.{^'n{\p+Vn? + \P-rin?) + ^^'n{\P+Vn? 

n 

- \V-^nf) + {iyn{p+rin){p-riny + C.C.)} (ic) 

F^ = Y.Xn\p^Vn\' (Id) 
n 

Here, ?7„, n=l,2,3, are the order parameter components, p± = {px ± ipy)/V^, Pa = —ida — 
2eAa/hc (e < 0), and a± = aoiT — T±) with T+ = Tcq + 2e/ao and T_ = Tco — e/oo- For 
the odd-in-T order parameter considered in the frequency domain, the rji are understood to 
be the leading order coefficients of the expansion about uj = We will restrict further 
analysis to the case e > (T+ > T_), which is a necessary condition to comply with the 
suggested T-breaking in phase B but not in phase A. The coefficients of the gradient terms 
are expanded to linear order in e and are defined by a number of phenomenological coupling 
constants Kj through = ki + K2e-Re(7"'), /x^ = K3e/m(7"'), i/„ = + K^e + ^667"", and 
Xn = 1^7 + Ks^R^il^) ■ The phase factor 7 = exp{i2TT/3) reflects the three-fold symmetry of 
the star of M. Note that the last term in Eq. (lb) favors T-breaking solutions for jS^ > 0. 

Model (1) accounts for a rich phase diagram already when H = 0. Fig. 3 summarizes 
the stability region of zero field solutions. Four regions in the parameter space can be 



distinguished according to their predicted phase sequence, as shown in Table 1 Most 
promising with regard to experiment are regions I and II. Both allow for two second order 
transitions at T+ and T^,(< T+), where the second phase manifestly breaks T-symmetry and 
can be identified with phase B of Fig. lb. 

For region I, we recover the experimental P — T diagram, when we assume a pressure 



dependent strain field, which decreases with increasing P and vanishes for P > P^ |2T[. On 
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the other hand, a third low-temperature phase for finite e is predicted for region II, which 
signals an unlocking of the relative phases of the order parameter components. Consequently, 
a fourth superconducting phase D appears between B and C, and all phase boundaries meet 
at a pentacritical point. Although this scenario seems to be at variance with Fig. lb, we 
cannot immediately reject it because the BD and DC boundaries are predicted to be very 
close together rendering an observation of the fourth phase difficult. Estimates for the model 
parameters, as derived from experimental data (shaded area inside the region of stability in 
Fig. 3) do not unambiguously favor one of the two regions. 

Let us now discuss the H — T diagram predicted by model (1). Generally, the upper 
critical field ifc2 is determined by the smallest eigenvalue of the linearized Ginzburg-Landau 



equations p!T| , p2| , which can be seen as a coupled system of Schrodinger-like equations, with 
the order parameter playing the role of a multicomponent wavefunction. In our case, this 
leads to three decoupled equations, which makes it possible to obtain all eigenvalues in a 
closed analytic form for a general magnetic field. For H = H[cosipsm.9,smipsm9,cosd), 
Hc2 is given by the largest of the three values 

Hn{T) = 

2|e| y/Rn + jJ,n COS 6^ 
Rn = (/in - knP) COS^ + Xnif^n ' i?e(z/„e*2^)) siu^ 9 (2) 

where = and ai = 02 = 

A tetracritical point in the H — T plane corresponds to a crossing of the two largest values 



Hn as a function of temperature. As mentioned by Garg p2[, such a crossing necessarily 
requires the existence of a conserved quantity associated with the Schrodiger equations. 
For our model, this property follows immediately from the decoupling of the equations 
dictated by translational invariance, therefore allowing, in principle, a tetracritical point 
for all directions of the magnetic field. Within a one domain model, the almost complete 



isotropy of Hc2 in the basal plane requires that the in-plane anisotropic gradient terms 
are very small {\vn\ ^ /^n)- For the extreme case Vn = 0, we find kinks in Hc2 for all field 
directions provided that ki > 0, K7 > 0, K2 + \k3\/^/3 > 0, and + Ki^g > p4|. These 



constraints do not represent strong limitations. Note, that the gradient terms proportional 
to e are necessary for a tetracritical point to appear because for vanishing strain the slopes 
of the Hn versus T lines would be identical. As a consequence, the slopes of the low-field 
and high- field i7c2-lines {AN and CN phase boundaries in Fig. la) approach the same value 
for e — 0, i.e. for P —* P^.. This holds for any field direction, and represents a qualitative 
difference from the predictions of the E model, where the slope difference approaches a 
nonzero limit for in-plane fields and vanishes for c-axis fields. 

Our model is further consistent with observed anisotropics of the penetration depth and 
Hcii and predicts a kink of Hd at for all field directions because of an onset of an 
additional order parameter component below T,,. 

The above given analysis of the phase diagram assumes a homogeneous symmetry break- 
ing field. This one-domain picture is an oversimplification of the situation in UPts. Any 
orthorhombic field is very likely to build a six-fold domain structure, as it was observed 
for the antiferromagnetic order with domain sizes of less than 150 A. Such an underlying 



domain structure can significantly alter the superconducting state |]25|. This point deserves 
further investigation to allow a more reliable comparison with experiment. 

Our model may resolve some additional phenomenological puzzles. Most notable is 
the absence of Josephson tunneling from conventional superconductors into heavy fermion 



materials, in particular for UPta [^]. Assuming the superconductivity to arise from FCM 
pairs with odd-in-T symmetry, no single pair Josephson coupling to a conventional even-in- 
T superconductor is possible without applied or internal (spontaneous) magnetic field that 
breaks T . Even in this case the staggered phase across an arbitrary sample face will yield 
zero net pair tunneling unless: (i) one cleaves along planes perpendicular to the pairing 
fc- vectors, or (ii) one uses atomic scale tunneling probes (e.g., STM), which can sample a 
small region of non-zero phase. This suggests a clear experimental test of the odd-in-T 
hypothesis. 
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grant from the Department of Energy, Office of Basic Energy Sciences, Division of Materials 
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TABLES 

TABLE I. Order parameters for stable phases of model (1) for H = depending on the 

parameter regions I-IV (cf. Fig. 3). Results apply for a constant (static) strain field e. Phase 1 
appears at T+, phase 2 at T* = r+ - - /32 + /33)(T+ - r_), < r+, and a possible third 

phase may exist below Tg = T+ - 2{/3i + p2 - p3)/{Pi - p2 - p3){T+ - r_), Tg < T*. x and ^ are 
temperature dependent real numbers, e = corresponds to the high-pressure phase (Fig. lb). 
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FIGURES 

FIG. 1. Schematic phase diagram of UPts in the P = and H = planes as deduced from 
experiments ||2|. Areas supposed to belong to the same phase in the three-dimensional H-P-T 
diagram are given the same letter. 

FIG. 2. kz = ^ cross-section of the Brillouin zone of UPts. kj are wave vectors of the star of 
the M-point representation used in our model. 

FIG. 3. Region of stability for solutions of model (1) for H = as a function of 
Ginzburg-Landau parameters for /?i > 0. The free energy has no lower bound in the unstable 
region, which holds also for (3i < 0. The shaded strip marks the region consistent with specific 



heat measurements [18,1^ (dotted lines: the specific heat jump ratio, r, at and with re- 
spect to the normal state, 1.25 < r < 1.33), estimates for T_ from //c2-measurements (dashed 
lines; x = (T^ — T_)/(T+ — T^), 0.2 < x < 0.45) and the assumption that only phase B violates 
T-invariance. Characteristic phase diagrams found in regions (I)- (IV) are described in the text 
and in Table 1. 
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Fig.2 




Fig.3 




